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Abstract 

We perform a detailed study of Gibbs-non-Gibbs transitions for the Curie- Weiss model 
subject to independent spin- flip dynamics ("infinite-temperature" dynamics). We show 
that, in this setup, the program outlined in van Enter, Fernandez, den Hollander and 
Redig 1^ can be fully completed, namely that Gibbs-non-Gibbs transitions are equivalent 
to bifurcations in the set of global minima of the large-deviation rate function for the 
trajectories of the magnetization conditioned on their endpoint. As a consequence, we 
show that the time-evolved model is non-Gibbs if and only if this set is not a singleton for 
some value of the final magnetization. A detailed description of the possible scenarios of 
bifurcation is given, leading to a full characterization of passages from Gibbs to non-Gibbs 
— and vice versa — with sharp transition times (under the dynamics Gibbsianness can be 
lost and can be recovered). 

Our analysis expands the work of Ermolaev and Kiilske who considered zero mag- 
netic field and finite-temperature spin-flip dynamics. We consider both zero and non-zero 
magnetic field but restricted to infinite-temperature spin-flip dynamics. Our results re- 
veal an interesting dependence on the interaction parameters, including the presence of 
forbidden regions for the optimal trajectories and the possible occurrence of overshoots 
and undershoots in the optimal trajectories. The numerical plots provided are obtained 
with the help of MATHEMATICA. 
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1 Introduction and main results 



Section 1.1 provides background and motivation, Section |1.2| a preview of the main results. 
Section |1.3| introduces the Curie- Weiss model and the key questions to be explored. Section lA 



recalls a few facts from large-deviation theory for trajectories of the magnetization in the 
Curie- Weiss model subjected to infinite-temperature spin-flip dynamics and provides the link 



with the specification kernel of the time-evolved measure when it is Gibbs. Section 1.5 states 
the main results and illustrates these results with numerical pictures. The pictures are made 
with MATHEMATICA, based on analytical expressions appearing in the text. Proofs are given 
in Sections |2| and ^. Section |^ takes up half of the paper. 

1.1 Background and motivation 

Dynamical Gibbs-non-Gibbs transitions represent a relatively novel and surprising phenome- 
non. The setup is simple: an initial Gibbsian state (e.g. a collection of interacting Ising spins) 
is subjected to a stochastic dynamics (e.g. a Glauber spin-flip dynamics) at a temperature 
that is different from that of the initial state. For many combinations of initial and dynamical 
temperature, the time-evolved state is observed to become non-Gibbs after a finite time. Such 
a state cannot be described by any absolutely summable Hamiltonian and therefore lacks a 
well-defined notion of temperature. 

The phenomenon was originally discovered by van Enter, Fernandez, den Hollander and 
Redig Q for heating dynamics, in which a low-temperature Ising model is subjected to 
an infinite-temperature dynamics (independent spin-fiips) or a high-temperature dynamics 
(weakly-dependent spin-flips). The state remains Gibbs for short times, but becomes non- 
Gibbs after a finite time. Remarkably, heating in this case does not lead to a succession 
of states with increasing temperature, but to states where the notion of temperature is lost 
altogether. Furthermore, it turned out that there is a difference depending on whether the 
initial Ising model has zero or non-zero magnetic field. In the former case, non-Gibbsianness 
once lost is never recovered, while in the latter case Gibbsianness is recovered at a later time. 

This initial work triggered a decade of developments that led to general results on Gibb- 
sianness for small times (Le Ny and Redig |l^, Dereudre and Roelly Q), loss and recovery 
of Gibbsianness for discrete spins (van Enter, Kiilske, Opoku and Ruszel |1C, ^ 



Redig, Roelly and Ruszel U^]), and loss and recovery of Gibbsianness for continuous spins 
(Kiilske and Redig |12], Van Enter and Ruszel (5|, ^). A particularly fruitful research direc- 



tion was initiated by Kiilske and Le Ny |^], who showed that Gibbs-non-Gibbs transitions 
can also be defined naturally for mean-field models, such as the Curie- Weiss model. Precise 
results are available for the latter, including sharpness of the transition times and an explicit 
characterization of the conditional magnetizations leading to non-Gibbsianness (Kiilske and 
Opoku [0], Ermolaev and Kiilske 0]). In particular, the work in 0] shows that in the mean- 
field setting Gibbs-non-Gibbs transitions occur for all initial temperatures below criticality, 
both for cooling dynamics and for heating dynamics. 

The ubiquitousness of the Gibbs-non-Gibbs phenomenon calls for a better understanding 
of its causes and consequences. Unfortunately, the mathematical approach used in most 
references is opaque on the intuitive level. Generically, non-Gibsianness is proved by looking 
at the evolving system at two times, the inital and the final time, and applying techniques from 
equilibrium statistical mechanics. This is an indirect approach that does not illuminate the 
relation between the Gibbs-non-Gibbs phenomenon and the dynamical effects responsible for 
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its occurrence. This unsatisfactory situation was addressed in Enter, Fernandez, den Hollander 
and Redig where possible dynamical mechanisms were proposed and a program was put 
forward to develop a theory of Gibbs-non-Gibbs transitions on purely dynamical grounds. 
The present paper shows that this program can be fully carried out for the Curie- Weiss model 
subject to an infinite-temperature dynamics. 

In the mean-field scenario, the key object is the time-evolved single-spin average condi- 
tional on the final empirical magnetization. Non-Gibbsianness corresponds to a discontinuous 
dependence of this average on the final magnetization. The discontinuity points are called bad 
magnetizations (see Definition below). Dynamically, such discontinuities are expected to 
arise whenever there is more than one possible trajectory compatible with the bad magnetiza- 
tion at the end. Indeed, this expectation is confirmed and exploited in the sequel. The actual 
conditional trajectories are those minimizing the large-deviation rate function on the space of 
trajectories of magnetizations. The time-evolved measure remains Gibbsian whenever there 
is a single minimizing trajectory for every final magnetization, in which case the specifica- 



tion kernel can be computed explicitly (see Proposition 1.4 below). In contrast, if there are 
multiple optimal trajectories, then the choice of trajectory can be decided by an infinitesimal 
perturbation of the final magnetization, and this is responsible for non-Gibbsianness. 

1.2 Preview of the main results 

In the present paper we study in detail the large-deviation rate function for the trajectory 
of the magnetization in the Curie- Weiss model with pair potential J > and magnetic field 
/i G M (see ( |1 . 1| ) below). We exploit the fact that, due to the mean-field character of the 
interaction, this rate function can be expressed as a function of the initial and the final 



magnetization only (see Proposition L2 below), i.e., the trajectories are uniquely determined 
by the magnetizations at the beginning and at the end (see Corollary |1.3| and Proposition 
below). Here is a summary of the main results (see Fig. ||): 

1. If < J < 1 (supercritical temperature), then the evolved state is Gibbs at all times. 
On the other hand, if J > 1 (subcritical temperature) there exists some time ^jj at 
which multiple trajectories appear. The associated non-Gibbsianness persists for all 
later times when h = (zero magnetic field). All these features were already shown by 
Ermolaev and Kiilske Q. 

2. For h ^ there is a time > at which Gibbsianness is restored for all later times. 

3. There is a change in behavior at J = |. For 1 < J < |: 

(a) If /i = 0, then only the zero magnetization is bad for t > '^c- 

(b) If h > (h < 0), then there is only one bad magnetization for ^jj < i < ^'^<. This 
bad magnetization changes with t but is always strictly negative (strictly positive). 

For J > |: 

(a) If h = 0, then there is a time > such that for < t < there are 
two non-zero bad magnetizations (equal in absolute value but with opposite signs), 
while for t > ^'c only the zero magnetization is bad. 

(b) If h ^ and small enough, then there are two times > between ^jj and 
such that for < t < and < t < 'f* only one bad magnetization occurs, 
while for < t < two bad magnetizations occur. 
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h = »— 



hj^O • • • • 

Figure 1: Crossover times for h = and h when J > |. 

All the crossover times depend on J, h and are strictly positive and finite. Our analysis 
gives a detailed picture of the optimal trajectories for different J, h and different conditional 
magnetizations. Among the novel features we mention: 

(1) Presence of forbidden regions that cannot be crossed by any optimal trajectory. The 
boundary of these regions is given by the multiple optimal trajectories when bifurcation 
sets in. The forbidden regions were predicted in and first found, for /i = 0, by 
Ermolaev and Kiilske 0. 

(2) Existence of overshoots and undershoots for optimal trajectories for h ^ 0. 

(3) Classification of the bad magnetizations leading to multiple optimal trajectories. These 
bad magnetizations depend on J, h and change with time. 



1.3 The model 

1.3.1 Hamiltonian and dynamics 

The Curie- Weiss model consists of N Ising spins, labelled i = 1, . . . , N with N gN. The spins 
interact through a mean-field Hamiltonian — that is, a Hamiltonian involving no geometry 
and no sense of neighborhood, in which each spin interacts equally with all other spins — . 
The Curie- Weiss Hamiltonian is 

N N 

H^{a) := aiaj -h^Ui, ae^N, (1.1) 

i,j=l i=l 

where J > is the (ferromagnetic) pair potential, /i G M is the (external) magnetic field, 
Otv := { — 1,+1}^ is the spin configuration space, and a := ((Tj)^]^ is the spin configuration. 
The Gibbs measure associated with is 

with the normalizing partition sum. 

We allow this model to evolve according to an independent spin-flip dynamics, that is, a 
dynamics defined by the generator L^r given by (see Liggett [0] for more background) 

N 

(L;v/)(^):= J^[/(0-/(a)], f-.n^^R, (1.3) 

1=1 
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where a* denotes the configuration obtained from a by fiipping the spin with label i. The 
resulting random variables a{t) := {(Ji{t))^i constitute a continuous-time Markov chain on 
J^TV- We write fi^ to denote the measure on Q]\i at time t when the initial measure is fi^ and 
abbreviate i^t '■= {fJ-t^)N€'N- 



1.3.2 Empirical magnetization 

To emphasize its mean- field character, it is convenient to write the Hamiltonian ( |1.1| ) in the 
form 

H^{a) = NH{mN{a)) (1.4) 

where 

H[x) := -\Jx'^ -hx, xG R. (1.5) 

and 

1 ^ 

mNi(T) := —^at (1.6) 

i=l 

is the empirical magnetization of cr G 17 jv, which takes values in the set A4n '■= {— Ij— 1 + 
2N~^ , . . . , +1 — 2A^^^,+1}. The Gibbs measure on induces a Gibbs measure on A^at 
given by 

/ N \ e~^-^("') 
~^^'^y-=[u^N)-^^' mGA^^, (1.7) 

where is the normalizing partition sum. 

The independent (infinite-temperature) dynamics has the simplifying feature of preserving 
the mean-field character of the model. In fact, the dynamics on O^v induces a dynamics on 
A^AT, which is a continuous-time Markov chain {m^)t>o with generator Lj\[ given by 

{LNf){m) := N[f{m - 2N-^) - f{m)] + N[f{m + IN'^) - f{m)] , (1.8) 

for /: A^AT — )• M. Adapting our previous notation we denote fi^ the measure on Af^r at 
time t, and abbreviate Jit '■= {fJ't^)NeN- Due to permutation invariance, characterizes 
and vice versa, for each and t. We write to denote the law of (m^)j>o, which lives on 
the space of cadlag trajectories L'[o,oo)([— 1) +1]) endowed with the Skorohod topology. 



1.3.3 Bad magnetizations 

Non-Gibbsianness shows up through discontinuities with respect to boundary conditions of 
finite- volume conditional probabilities. For the Curie- Weiss model it is enough to consider 
the single-spin conditional probabilities 

7^ (-71 I aN-i) ■■= /if (<7l I (JN-l) , (1.9) 
defined for ai G { — 1, +1} and a]\f~i G AIat-i, and any spin configuration (Tiv_i G ^n-i such 



that uiN-iicTN-i) = OiN-i- By permutation invariance, ( |l.9[) does not depend on the choice 
of aN-i- 

The central definition for our purposes is the following. 



5 



Definition 1.1. (Kiilske and Le Ny g) Fix t > 0. 

(a) A magnetization a S [— 1,+1] is said to he good for nt if there exists a neighborhood Ma 
of a such that 

jt{- I a) := lim 7^ (• I a^-i), (1-10) 

exists for all a in J\fa and all {aN)NeN such that aN & -Mn for all N G N and liniAr^oo ctN = 
a, and is independent of the choice of {a]\f)NeN- The limit is called the specification kernel. 
In particular, a 1— )• 7t(- | a) is continuous at a = a. 

(b) A magnetization a £ [— 1,+1] is called bad if it is not good. 

(c) fit is called Gibbs if it has no bad magnetizations. 

1.4 Path large deviations and link to specification kernel 

The main point of our work is our relation between path large deviations and non-Gibbsianness. 
For the convenience of the reader, let us recall some basic large deviation results for the Curie- 
Weiss model. For background on large deviation theory, see e.g. den Hollander [^. 

1.4.1 Path large deviation principle 

Let us recall that a family of measures on a Borel measure space satisfies a large deviation 
principle with rate function / and speed if the following two conditions are satisfied: 

liminf — logi^^(^) > — inf I(x) for A open (l-H) 
lim sup log i/^ (^4) < — sup/(x) foreclosed (1-12) 

The proof of the following proposition is elementary and can be found in many references. 
The indices S and D stand for static and dynamic. 

Proposition 1.2. (Ermolaev and Kiilske 0, Enter, Fernandez, den Hollander and Redig |^) 

(i) (/i^)iVGN satisfies the large deviation principle on [— 1,+1] with rate N and rate function 
Is — inf(/5') given by 

— — — 1-1- m 1 — m 
Is{m) := H{m) + I{m), I{m) := — ^ — log(l m) H ^ — log(l - m). (1.13) 

(a) For every T > 0, the restriction of {P^)]sf^^ to the time interval [0,T] satisfies the large 
deviation principle on D^Q^rp^{[—l, +1]) with rate N and rate function I^ — inf(/'^) given by 

I^{ip):=Is{m)+ll{^). (1-14) 

where 

jT^^._(jQL{ip{s),ip{s))ds if (f exists, 
^ \ 00 otherwise, 

is the action integral with Lagrangian 



1 r-, o 1 , / \/4 (1 - m2) + rh'^ + m\ 

L{m,m) := --V4(l - m?) + Tf)T? + -mlog ( ^ — _ 1 + 1. (1.16) 
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Let 

g^^(m) := (m7v(0) = m I m7v(t) = a) , m £ Mn (1-17) 

be the conditional distribution of the magnetization at time given that the magnetization 
at time t is a. The contraction principle applied to Proposition |1 .21 (11) implies the following 
large deviation principle. 

Corollary 1.3. For every t > and a G [—1,+!], {Qt^a)NeN satisfies the large deviation 
principle on [—1,+!] with rate N and rate function Ct,a — inf(Ct^a) given by 

Ct,a{m):= inf l\ip). (1.18) 

ip{t)=a 

Note that 

inf Cta{m)= inf inf /*((/?)= inf P{ip). (1.19) 

mG[-l,+l] ' »ne[-l,+l] tp: ip(a)=m, p: 'p{t)=a 

ip{t)=a 

1.4.2 Link to specification Icernel 

The following proposition provides the fundamental link between the specification kernel in 
(pop and the mimmizer of ( p^ when it is unique, and is a straightforward generalization 



to arbitrary magnetic field of a result for zero magnetic field stated and proved in Ermolaev 
and Kiilske ffTtl. 



Proposition 1.4. Fix t > and a G [— 1,+1]. Suppose that ( |1.19 ) has a unique minimizing 
path {ipt,ais))o<s<t- Then the specification kernel equals 

7tiz \ a) = fy \ , ZG{-1,+1}, (1.20) 

where pt{-,-) is the transition kernel of the continuous-time Markov chain on { — 1,+1} jump- 
ing at rate 1, given by pt{l,l) = pt{—l,—l) = e~* cosh(t) and pt{—l,+l) = pt{l,—l) = 
e~* sinh(t). 



Remark: Note that the expression in the right-hand side of (1.20) depends on the optimal 



trajectory only via its initial value (^1,0 (0). Thus, ( |L20|) has the form 



-ft{z\a) = rt{z,J^t,a{0) + h), (1.21) 

where (pt,a{0) is the unique global minimizer of m 1— t- Ct^Q,(m) and m 1— t- Tt{z, m) is continuous 
and strictly increasing (strictly decreasing) for 2; = 1 (2 = — 1). 

1.4.3 Reduction 

The next proposition allows us to reduce ( |1.19| ) to a one-dimensional variational problem. 
Consider the equation 

k-^^^{m) = lt,a{m) (1.22) 

with 

k-^^^{m) := a\m) cosh(2/i) + b-^{m) sinh(2/i), 
k,a{n^) ■= mcoth(2t) — acsch(2t). 
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where 

a"^ (m) := sinh(2 Jm) — m cosh(2 Jm), 
{m) := cosh(2 Jm) — m sinh(2 Jm). 

Proposition 1.5. Let Ct^a be as in ( |1.18| ). Then, for every t>0 and a G [— 1,+1], 



Ct,a{m) = Is{m) 
1 

+ - 



+2 



4t + log 
a log 



1 2 + log 

1 — m^y 
- Cie-^^ + Cse^* 
1 — a 



1-R- 2Ciae-2* 
1 + - 2C7iae-2* 



1 + i? - 2Cim' 
1-R- 2Cim 



mloG 



- Ci + C2 
1 — m 



C2 



Ci{t, a, m) 
C2(t, a, m) 

ii(Cl,C2) 



Q— me 



= V1-4C1C2. 

Furthermore, the critical points of Ct,a ol"^^ the solutions of ( |1.22 ) . Hence, 



inf /*((/?) 



mirir 



Ct,a{m) , 



ip: ip{t)=a ' ' TO solves (1.22) 

and the constrained minimizing trajectories are of the form 

ipf^^{s) := csch(2t)|msinh(2(t - s)) + asinh(2s)| 



<s<t 



rh = rh{t, a) 
The identities 



argmiii 



solutions of ( 1.22 ) 



and 



k-^'^{m) = 2 cosh^ (Jm + h) [tanh( Jm + h) - m] + 



lim It o(m) = m 

t— >-oo ' 



m 



(1.24) 



(1.25) 



(1.26) 



(1.27) 

(1.28) 
(1.29) 

(1.30) 
(1.31) 



imply that in the limit t — t- cxd ( |1.22D reduces to tanh( Jm + h) = m. This is the equation for 
the spontaneous magnetization of the Curie- Weiss model with parameters J, h. This equation 
has always at least one solution and the value 

m°° = m^{J, h) := the largest solution of the equation tanh( Jm + h) = m (1.32) 

is well known to be strictly positive if /i > or if J > 1. In these regimes, the standard 
Curie- Weiss graphical argument shows that, for m > 0, 



m = m 



(1.33) 



We also remark that when t — t- the function It^a converges to the line defined by the equation 
m = a. This implies that for short times there is a unique solution of ( |l.22D and it is close to 
a. 
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1.5 Main results 



In Section |1.5.1| we state the equivalence of non-Gibbs and bifurcation that lies at the heart 
of the program outlined in (Theorem 1.6). In Section 1.5.2 we introduce some notation. 
In Section 1.5.3 we identify the optimal trajectories for a = 0, h = (Theorems 1.7-L^). 
In Section 1.5.4 we extend this identification to a S [— 1,+1], /i € M (Theorem |1.9|). In 
Section 1.5.5 we summarize the consequences for Gibbs versus non-Gibbs (Corollary |1.1C1|). 



1.5.1 Equivalence of non-Gibbs and bifurcation 

The following theorem proves the long suspected equivalence between dynamical non-Gibbsianness, 
i.e., discontinuity of a i— t- jti' \ ct) at qq, and non-uniqueness of the global minimizer of 
m I—)- Cf^Q,Q(m), i.e., the occurrence of more than one possible history for the same a. 

Theorem 1.6. a i— )• 7t(o" | a) is continuous at oq if and only i/inf^. ip(^t)=ao iv^) ^'^^ 
unique minimizing path or, equivalently, inf^gj.^ Ct^a(j(m) has a unique minimizing mag- 
netization. 



1.5.2 Notation 



Due to relation ( |1.27 ), our analysis focusses on the different solutions of ( [1.22 ) obtained as 



a are varied. In particular, we must determine which of them are minima of the variational 
problem in ( |1.19| ). We write 

^t,a '■= the set of global minimizers of Ct^a- (1-34) 

For brevity, when a is kept fixed and At^a is a singleton {m(t,a)} for each t, we write rh{t) 
instead of m{t,a). When /i,a = 0, by symmetry we have A^^o = {0} or A^^^ = {i^(0}i 
where in the last case we denote by fh{t) the unique positive global minimizer. If both the 
initial and final magnetizations are fixed, then there is a unique minimizer that we denote as 



in ( |L28|) . That is, 

^T,^:= argmin I'M (1.35) 

if: ip{Q)=m, 
ip{t)=a 

for m,a £ [— 1,+1]. We emphasize that, by definition, Ct^ai^) = I^iv'Ta) 0t,ais) = 
(p^^'°'\s),s G [0, t]. In particular m{t,a) = (pt,a{^)- 

1.5.3 Optimal trajectories for a = 0, /i = 

The following theorem refers to a critical time 



i arccoth(2 J - 1) if 1< J < ^ 
if J > 3 



m, = m,{J):={^ ' 3"" (1-36) 

2' 



where t* = t*(J) is implicitly calculable: = t{'m^) where the function t{m) is defined in 
(pTl| ) below and rritf = rritf^J) is the solution of (|2.1^ ). 



Theorem 1.7. (See Fig. 0.) Consider a = and h = 0. 



9 



(i) IfO<J< 1, then 



(a) // 1 < J < |, then 



At,o = {0}, Vt > 0. 



{0} tfO<t<^„ 
{±m{t)} ^ft>^c, 



where 1 1— )• m{t) is continuous and strictly increasing on [^c, oo) with m(^c 



(Hi) If J > |, then 



At 



{0} 



ifO<t<^c, 

'c, 



{±m(t)} ift>^c 
where 1 1— )■ rh{t) is continuous and strictly increasing on [4'c,oo) with 7h{^c) 



(1.37) 
(1.38) 

= 0. 

(1.39) 
= : > 0. 



m(t,aj 



m(t.a) 



m(t,a) 



t, 'P, 




t{= tl) < 



t = ^r 



ti= t2) > ^-c 



Figure 2: Illustration of Theorem |l.7| . First row: Time evolution of the minimizing trajectories 
±{iptfl{s))o<s<t for t < ^'c, t = and t > for an initial Curie- Weiss model with (J, h) = 
(1.6,0) [regime (iii) in the Theorem]. The shaded cone is the forbidden region. Second row: 
Plot of m I— )• Ctfi{m) for the same times and parameter values. 



Let At^o('^) denote the cone between the trajectories ±ipt,Q- As a consequence of the 
previous theorem, no minimal trajectory conditioned in t' with t' > t can intersect the interior 
of this region. Such a cone corresponds, therefore, to a forbidden region. Forbidden regions 
grow, in a nested fashion, as the conditioning time t grows. There is, however, a distinctive 
difference between regimes (ii) and (iii) in the previous theorem: In Regime (ii) the forbidden 
region opens up continuously after while in Regime (iii) it opens up dis continuously. These 
facts are summarized in the following theorem. 

Theorem 1.8. Suppose that a = and h = 0. 

(i) J I— 7- m^,{J) is strictly increasing on (|,oo). 

(ii) J I— )• ^dJ) is strictly decreasing on (l,oo). 

(iii) J I— 7- Atfi{J) is left- continuous at J = ^ for all t > ^c(|). 

(iv) J I—)- A^^(-j) o(^) is right- continuous at J = J for all J > |. 

(v) For every J < 3/2 the map 1 1— )• Atfi{J) is continuous. 
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(vi) For every J > 3/2 the map t i— )■ Atfi{J) is continuous except at t = '^c. where it exhibits 
a right- continuous jump. 



1.5.4 Optimal trajectories for a G [— 1,+1], ^ G M 
Fore fixed (J, h) and a we say that there is (See Fig. ^): 

• No bifurcation if A^^q = {m(t,a)}, for all t > and the map 1 1— >• m{t,a) is continuous 
on [0, oo). 

• Bifurcation when there exists a < < oo such that 1 1-)- m{t, a) continuous except at 
t = tB and \Atg^a\ = 2. 

• Double bifurcation if there exist times < sb < ts < oo such that 1 1— t- m(t, a) continu- 
ous except at t = tB and t = sb, and |As^^q,| = |Atg^Q,| = 2. 

• Trifurcation if there exists a < < co such that t i— )• a) is continuous except at 
t = tx and |Atj,^Q,| = 3. 

The bifurcation times tB and sb, the trifurcation time and the trifurcation magnetization 
Mt (defined below) all depend on J, h. 



a=0.057,?a=0.0600321 



a=-OA 43,^B=0.259404,.?„=0.072933 1 
C 



Q'=-0.0761756,?a=0.0886908 




Bifurcation 



Double bifurcation 



Trifurcation 



Figure 3: Different scenarios for the evolution in time of m i— )• Ct^ai^)- Drawn lines: t = tB, 
t = tB, SB, t = tx (times at which multiple global minima occur or, equivalently, discontinuity 
points of 1 1— )• m(t,a)). Dotted lines: earlier time. Dashed lines: later time. 

The following theorem summarizes the behaviour of At^a (and therefore of the minimizing 
trajectories (ft, a) for different t,a. For J > |, let 



F{m) :-- 



csch[arccoth(A;'"''^(m))] 



Ub = UB{J,h) := max F{m), 

me[0,l] 

Lb = Lb{J, h) := min F{m). 

me[-l,0] 



(1.40) 



Theorem 1.9. (See Figs. 

(1) Suppose that y^^{a) / 0. 

(la) If k^'^{a) > and a > 0, then there are > and tR = tji{m^) > (implicitly 
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calculable from ( |3.8D ) such that t i— )■ m(t) is strictly increasing on [0, tf>] and strictly 
decreasing on [tji, oo) with rh{tji) = > m°°. 

(lb) If k^'^{a) < and a > 0, then 1 1— )• m(t) is strictly decreasing on [0, oo). 

(Ic) If k'^'^(a) > and a <0, then 1 1— )• rh{t) is strictly increasing on [0,oo). 

(Id) If k'^'^{a) < and a < 0, then there are > and tji = tji{m^) > (implicitly 

calculable from ( p.9D ) such that t i— )• rh{t) is strictly decreasing on [0, t/j] and strictly 

increasing on [tji,oo) with rh{tii) = fnfi < a. 

In all cases m(0) = a and linit^oo rh{t) = m°° . 

(2) Suppose that h = 0. 

(2a) 7/0 < J < 1, then there is no bifurcation. 

(2b) 7/ 1 < J < |, then there is bifurcation only for a = 0. 

(2c) If J > |, then there is bifurcation if a ^ {—Ub,Ub) and no bifurcation otherwise. 

(3) Suppose that h > 0. 

(3a) 7/0 < J < 1, then there is no bifurcation. 

(3b) 7f 1 < J < |, then there is bifurcation for a G [—1,Ub) and no bifurcation for 
a G [Ub,1]- 

(3c) If J > ^, then there exists a h^, = h^,{J) > such that 

- for every < h < h^, there exists a Mt G {Lb,Ub) with Mt < such that there is 

* no bifurcation for a G [iJe , 1] , 

* bifurcation for a G (Mp, Ub), 

* trifurcation for a = Mt, 

* double bifurcation for a G {Lb, Mt), 

* bifurcation for a G [— 1, L^]. 

- for every h > h* the behavior is the same as in (3b). 

In all cases a i— )■ t_B(a) is continuous and decreasing and a i— )• ^^(a) is continuous and 
increasing. 



Theorem 1.9 gives a complete picture of the bifurcation scenario. Regime (1) — which includes 
cases with zero and nonzero magnetic field — describes two types of behavior of optimal 
magnetization trajectories: monotone trajectories [cases (lb) and (Ic)] and trajectories with 
overshoot [cases (la) and (Id)]. In the latter, rh{t) increases to some magnetization larger 
(m^ smaller) than m°° and afterwards decreases (increases) to m°°. Regimes (2) and (3) refer 
to the existence of bifurcations and trifurcations. We observe that the different bifurcation 
behaviors — no bifurcation, single and double bifurcation — hold for whole intervals of the 
conditioning magnetization. In contrast, trifurcation appears at a single final magnetization 
for each h ^ 0. 

1.5.5 Gibbs versus non-Gibbs 



Theorem |1.6| establishes the equivalence of bifurcation and discontinuity of specifications, as 
proposed in the program put forward in [|3[. Due to this equivalence, the following corollary 
provides a full characterization of the different Gibbs-nonGibbs scenarios appearing during 
the infinite-temperature evolution of the Curie- Weiss model. Let 

< ^[7 := tBiUB) < := tB{MT) < ■■= tB{LB) < ^* := tB{-l), (1.41) 

and let Mb be the solution of tB^MB) = ^l- Denote Vt C [— 1,+1] the set of a-values for 
which a !->• 7t(-|a) is discontinuous. 
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m(t,a) 

1.0 r 



rh(t,a) 

1.0 r 



0.5 / 1.0 



Regime (la) 



m(t,a) 

1.0 r 



Regime (lb) 



rh(t,a) 

1.0 r 



a 



Regime (Id) 



a 



a 



Regime (Ic) 



Figure 4: Different regimes of Tlieorem |1.9| . Evolution in time of the minimizing trajectories 
M'Pt,a{s))Q<s<t for t <tR (dotted), t = Ir (drawn), t > tn (dashed). 
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Corollary 1.10. (See Fig. |.) 

(1) Leth = 0. 

(la) //O < J < 1, the evolved measure fit is Gibbs for all t >0. 
(lb) // 1 < J < |, then fit is 

- Gibbs forO<t< 

- non-Gibbs for t > with Vt = {0}. 
(Ic) If J > ^, then fit is 

- Gibbs forO<t< ^u, 

- non-Gibbs for t > ^jj with 

* Vt = {±a] for some a G {-Ub, Ub) if-^u <t< 
*A = {0} ift>^,. 

(2) Let h>0. 

(2a) // < J < 1, then fit is Gibbs for t > 0. 
(2b) // 1 < J < |, then fit is 

- Gibbs forO<t< "^u, 

- non-Gibbs for "iffj < t <"i!^ with Vt = {a} for some a G [—1, Ub), 

- Gibbs for t > 

(2c) If J > ^ and h < h* small enough, then fit is 

- Gibbs forO<t< "^u, 

- non-Gibbs for ^jj < t < with 

* T>t = {a} for some a £ [Mb, Ub) if^u<t<^L, 

* T>t = {a^,a'^} for some a^a^ G {Lb, Mb) if^L<t< ^t, 

* Vt = {a} for some a G [-1,Mt] if'^T<t<^*. 

- Gibbs for t > If h > h* , then the behaviour is as in (2b). 

In all cases Q^,c?,a depend on {t,J,h). 

2 Proof of Proposition |1.5| and Theorems |TTB|-|TT^ 



Proposition 1.5 is proven in Section 2.1, Theorems 1.6-1.8 are proven in Sections 2.2-2.4 



2.1 Proof of Proposition |1.5 

Proof. First note that, by ( 1.14| ), 



inf l\(f) 

ip: 'p{t)=a 



inf 

mGf-l,+ll 



Is{m)+ inf Ijjiif) 

tp: ip(0)=in, 
<p(t)=a 



> . 



(2.1) 



It follows from ( [L.14 - [lT5 ) and the calculus of variations that the stationary points of the 
right-hand side of (|2.lD are given by the Euler-Lagrange equation, complemented with a free- 
left-end condition and a fixed-right-end condition: 

— \ 
dm 

s=o dm 
<f{t) = a. 



^|^(v5(s),'/j(s)) 
as am 

^{^{s),^{s)) 



-{ip{s),Lp{s)), 



sG (0,t), 



(2.2) 
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Figure 5: Summary of Corollary l.lCj : Time versus bad magnetizations for different regimes. 
On the vertical a-axis, indicated by a thick line, is the set of bad magnetizations. G=Gibbs, 
NG=non-Gibbs. 



The first and the third equation i n (p.2| ) come from the third infimum in (2.1) and, together 
with (1.16), determine the form ( [1.28[ ) of the stationary trajectory. Inserting this form into 
( 1.14 ) we identify 

^*(Ca) = CtArri), (2.3) 

as stated in (|1.25| )~( 1.26 ). This identity reduces ( 1.19 ) to a one-dimensional variational prob- 
lem, 

inf I\ip)= inf I*(CJ= inf Ci,„(m) (2.4) 

(fi: ^{t)=a mG[-l,+l] me[~l,+l] 



The second equation in ( |2.2D corresponds to the second infimum in ( |2.1[) or, equivalently, to 
the rightmost infima in ( p.4[ ). It gives a trade-off between the static and the dynamic cost, 
establishing a relation between the initial magnetization and the initial derivative. After some 
manipulations this equation can be written in the form 



\q = a'\m)cosh.{2h) +h'\m)siTih{2h), m = (^^;*„(0), g = (^™„(0). (2.5) 



Differentiating ( 1.28 ), we get 



"ft 



*^(s) = 2csch(2t)|acosh(2s) - mcosh(2(t - s))|. 



(2.6) 
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and eliminating q from the last identity in (^]^) in favor of t and a, we conclude that m must 
be a solution of ( 1.22| ). Imposing this restriction to the chain of identities ( |2.4D we obtain 



([L2^ ) and hence (|L2^)-([T:29D . □ 



Prom now on our arguments rely on the study of ( 1.22| ), combined with continuity prop- 
erties of Cj Q as a function of t, a. 



2.2 Proof of Theorem 1.6 



Equation ( |l.20D follows in the same way as in the proof of Theorem 2.5 in Having disposed 
of this identity, we can now proceed to prove the equivalence. The proof relies on the following 
lemma. 

Lemma 2.1. For any t > and G [— 1,+1], there exists an open neighbourhood Na^ 7^ 
of the later, such that for all a G AAqq \ {oq} 

1. Ct^a has only one global minimum, namely, rri{t,a). 

2. a ^ rh{t,a) is continuous at a. If Ct^ao has a unique global minimum, the continuity 
is also valid at a = uq. 

3. If Ct^ao has multiple global minima, for two of them, namely, 771^(^,00) and mB{t,ao) 

lira rh{t, a) = rfiAit, uq) and lira 7h{t, a) = Tfisit, uq). 



Proof. A straightforward study of ( [1.22 ) shows that Ct^a has a finite number of critical points, 
for every fixed choice of J,h,t,a. 

Clearly, a 1— )• Ct^a and a 1— )• It^a are continuous with respect to the infinity norm in 
C([— 1, +1],]R). This, together with the fact that the left-hand side of ( 1.22|) does not depend 
on a, implies continuity of any critical point with respect to a. 

Let mi{t,ao), i = 1,. . . ,v be the global minima of Ct,ao- By continuity of the critical 
points, there exists a neighbourhood Mao ^"^^ smooth functions Mao ^ Oi ^ rni{t,a), i = 
1, . . . ,v, such that 

i. mi{t,a) are local minima of Ct,a, 
a. lim rni(t,a) = 7ni{t,ao). 
This properties proves the lemma if v = 1. Otherwise, let 

Bi{a) := Ct,a{rni{t,a)). 

The minimal cost is attained at the smallest of them: 

Ct,a{rh{'t,Oi)) = raiaBi{a) . 

Note that there is coincidence at oq due to the assumed multiplicity of minima: 

B{ao) := Bi{ao) , i = l,...v. (2.7) 
We expand the functions Bi up to first order order 

Bi{a) = B{ao) + S-(ao)(a - oq) + 0{a - oq) , (2.8) 
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and observe that, 

Biiao) ^ B'jiao), i^j. (2.9) 
The latter is due to the strict monotonicity of and the fact that each a) is a critical 



points of the function Ct^a{ ■ )• From (|2.2| )~ (|2.9| ) we conclude that for a in a possibly smaller 
neighbourhood ^ Mao there is a unique global minimum, and that property 3. holds with 

a = argminSj'(ao)) b = argmaxi?j'(ao), 
i i 

and 

mA{t, ao) := rha{t, ao), rnB{t, ao) := rnh{t, uq). 

□ 



We are ready to prove Theorem 1.6. 



Proof. Suppose that Ct^ao has a unique minimizer, denoted by ao) and let be the 
neighbourhood of the previous lemma. Then ( |1.21| ) holds for every a E Mao -> the continuity 
of m I— )• Tt{z, m) for every t, z gives the desired continuity of a i— )• 7f(- | a) at a = oq. 

To prove necessity, assume that Ct^ao has multiple global minima. Consider tjia and ms 
as in the previous lemma. Then, we have that there exist sequences < oq < converging 
to ao and such that 7t(- | a^) = Tt{-,m{t,a^)) and 

lim ■m{t,a~) = mB{t,ao) mA{t,ao) = lim m{t,a^). (2.10) 

n— >oo n^oo 

Again using continuity of Tt with respect to m, we get 

lim 7t(z I a~) = Tt{z,mB{t,ao)) ^ Tt{z,rhA{t, uq)) = lim 74(2: | a+). 

n— >oo n— >oo 

Hence ao is a bad magnetization. □ 



2.3 Proof of Theorem 1.7 



To determine which solutions of (|1.22 ) are global minima of Ctfi when h = 0, we will pursue 



the following strategy. Using ( 1.22| ) we can write t as a function of 



m: 



t{m) := \ arccoth ( ^ J • (2.11) 

This allows us to determine for which time t the magnetization m can be a possible minimum 
(i.e., a solution of ( 1.22D ). 



Lemma 2.2. Let A C [— 1,+1] be the set of m -values such that m is the solution of (1-22) 
for some t > 0, i.e., A = {m G [—1, +1] : a'' {m)/m > 1}. Then, for every m & A, 

Ci(m),o("T') = ^Jm?' + i log [1 — mtanh(Jm)] =: CM{fn). (2.12) 



In words, (2.12) is the cost for m at the time at which it is a possible minimum. 



Proof Insert ( ^HD into ([L2^ ) and use (^). □ 
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We now start the proof of Theorem 



1.7. 



Proof. First note that Itfi is Hnear with slope coth(2t) G (l,oo) and lt,o{0) = 0, and that a'^ is 
antisymmetric. Hence, if m is a solution of ( |1.22| ), then also — m is a solution. Further note 
that 

a'^'(m) = (2 J — 1) cosh(2Jm) — 2Jmsinh(2Jm), 
a-^" [m) = 4:J{J - 1) sinh(2Jm) - 4J^m cosh(2 Jm). 



(i) If < J < ^, then a"^' {m) < for all m, and hence m = is the unique solution for all 
t > 0. If ^ < J < 1, then a'^" [m) < for all m, hence a'^ is convex, and so it suffices to 
compare slopes at 0: A;'^'°'(0) = a-^'(O) = 2 J - 1 < 1 and /t,o'(0) = coth(2t) > 1. Again, m = 
is the unique solution for alH > (see Fig. 

(ii) As before, a'^" {m) < for all m, but now the slopes at can be equal, which occurs 
when t = with *c defined in (|1.36| ). This proves that At = {0} for < t < and 
At C {— m(t), 0, = the set of solutions of ( |1.22| ) for t > ^I'c. It is easily seen from 
Fig. |7(a)| that rh{t) is continuous and strictly increasing on [^cj oo) and m(^c) = 0. It 
remains to show that {— m(t),m(t)} are the global minima for all t > ^c- This follows from 
the strategy behind the proof of Lemma 2^. Since Ctfl{0) = for all i > 0, it suffices to prove 
that m I—)- CMi'fn) is strictly decreasing. From (2.12) we have 



da 



M 



t(m),0 



dm 



m) + 



da 



t{m),0 



dt 



(m) t'{m). 



(2.14) 



The first term is zero by the definition of t(m) (each m is a stationary point of Ctfi at time 
t = t{m)). The second term is < because t'{m) > and 



dCt^ 
dt 



(m) = L(^-oW,Co(*)) 

t 

+ 







dt 



dm 



dt 



ds 



t 



dt 



L{0,r 



dL 



(0,r)r 



dm 

= -i\/4 + r2 + l<0 

(2.15) 

with r = ipY\){t), where the second equality uses ( p. 21 ). Since Cm(0) = 0, this yields the claim 
(see Figs. fTjaimMcli) ■ 

(iii) This case is more difficult, because a'^ no longer is convex on (0, 1). Let ^'J be the first 
time at which a solution different from exists. To identify ^'J, let 



Tm{x) := (x — m)a"^ (m) + a'^ (m), 



(2.16) 
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and let mi be the solution of the equation Tmi(O) = = — mia"^'(mi) + a"^(mi), i.e., 

a"^(?7ii) 



mi 



a-^'(mi) 



(2.17) 



From mi we get by using ( ^.11 ): = t{mi). As before, a solution of ( 1.22| ) for t > is 
not necessarily a minimum. To find out when it is, we follow the same strategy as in case (ii). 
Again, Ctfl{0) = for all t > 0, and hence we must look for > such that 



CM{m^) = 0. 

Knowing m^, , we are able to compute ^'c using ( p. Ill) , 

= t(m*). 



(2.18) 



(2.19) 



In words, is the first time at which no longer is a minimum. As in case (ii), it suffices to 
prove that m i— )■ Cm(?7i) is strictly decreasing on (m=K,oo). Again, we have (|2.14D . Since 



,, , 1, , /a'^("^)\ f J', X a-'fm)! 1 

*'(,„) = -(axctanh)' M-i (m) - ^\ -. 



(2.20) 



it follows that t'{m) = if and only of m = a"^ {m) / a"^ ' {m) , which is the same condition as 
( |2.17| ). This gives us a graphical argument to conclude that t'{m) < for < m < mi and 
t'(m) > for m > mi (see Figs. 7(d), 7(e), 7(f)| ). On the other hand, m^, > mi. □ 



1.0 

J=0.4 0.5 




-1.0 \ 


-0.5 


-0.5 
-1.0 


' ' ^ — \o 



Figure 6: m i->- {m), m i— )■ lt,o{m) for Regime (i). 
2.4 Proof of Theorem ITTS 



Proof. From ( |1.28| ) it follows that 
m > m, i > t 

(i) First note that, because Cm("^*) = 0, 

9m* 



Co(«)>Co(^) yo<s<t. 



dJ 



(m* 



M 



dm 
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(2.21) 



(2.22) 



CJm) CJm) 




Regime (iii), J = 2. 



Figure 7: (a+d) m i— )• a-^{m), m i— )• ltfl{m); (b+e) m i— )• Ctfl{m), for < t < (dotted), 
t = ^c (drawn), t > (dashed); (c+f) m i-)- CM{m). 



As in Section 2^, case (iii), we have {dC m I dm){m^) < 0. But 

dCM 



dJ 



-(m*) > 



< tanh(Jm*), 



(2.23) 



which yields the claim because < ?tt-°°. 

(ii) The claim is straightforward for 1 < J < |. For J > | we need to prove that the function 



J — )• a"^ {m^,)/m^ is strictly increasing. In fact, 



d 




1 




dJ 




ml 


dJ \ 



dm 



a7 



(2.24) 



= 1 + cosh(2Jm*) — sinh(2Jm*) . 

The strict positivity of the last expression is equivalent to the inequality m^, < coth(Jm*), 
which is satisfied for J > 1. 

3 3 

(iii) We have a (m) J, a 2 (m) as J J, | for all m G (0, 1), with a and a 2 continuous. By Dini's 
theorem, the convergence is uniform. 

(iv) Since ^^dJ) < ^c(|)i the same argument as in case (iii) can be used. 



(v) and (vi) are consequences of parts (ii) and (iii) of Theorem 1.7 



□ 



3 Proof of Theorem 

In Sections |3.1| and we prove that overshoots, respectively, bifurcations, take place in 



regime (1), respectively (2)-(3), of Theorem L£. The analysis of the former regime does not 



distinguish on whether the initial field h is zero or not. 
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3.1 Regime (1): Overshoots 

The trick is again to write f as a function of m. From ( |1.22[) , we have 

- k-^'^{m) + 2mcoth(2t) = mcoth(2t) +Qcsch(2t). 

Hence, from ( 1.23| ), 

t^'^{m)[-k-^'''{m) + 2mcoth(2t)] 

= [mcoth(2t) - acsch(2t)] [mcoth(2t) + acsch(2t)] 



which imphes 
Solving for t we find 



a 



[m 



Q^) coth^(2t) - 2mt^'^{m) coth(2t) 



(3.1) 



(3.2) 



(3.3) 



tpim) 



1 fmk ' {m) + \a\ J <^{jn) 
- arccoth 

2 V — 



1 
2 




arccoth 



\k'^'^(ra)\ ^ + 
2m k-^'^{m) 



arccoth 



— arccoth 
2 



/mk ' {m) — |q| y¥(m)\ 
V m? — o? / 



2m k-^'^{m) 



if m ^ iba 

if m = ±a , m k'^'^{m) < 
if m = iba , m k"''^{m) > 

if m 7^ iba 

if m = ±a , mk'^'^{m) > 
if m = iba , mk"^'^ (m) < 



(3.4) 



(3.5) 



with 



$(m) := k-^'^{mY - + 



(3.6) 



These are times at which m is a stationary point (not necessarily a minimum) both for Ct^ai'm) 
and for Ct-a{Tnn) [Equation (3^), and hence the solutions (^jj) are insensitive to the sign of 
a]. Overshoots and undershoots occur for values of m satisfying (i) tp{m) > and tL{m) > 
and (ii) at these times m is a minimum. 

The m-dependence oitp and t/, is depicted in Figs ^ and ^ for cases in which both functions are 
injective, i.e., for a for which there is only one critical point for each t. In more complicated 
cases, for instance, when overshoot and bifurcation occur simultaneously, there are two or 
more stationary points only one of which is a minimum. 

We divide the analysis in four steps. 

Step A: Existence of , and z ^. We observe that there exists a unique m°° > such 
that k'^'^{m°°) = m°°. Furthermore, 



k'^^'^im) = 



tanh(2/i) 



—a'^{m) 
{m) 



--: A{m) 



(3.7) 



This function A has the features depicted in Figure IC: it is odd, satisfies A{0) = 0, A{1) = 1 
and it is convex with only one global minimum between and 1. We conclude that there 
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Plot of tp (dotted), tL (dashed) for a and —a. Plot of i i-^- rh{t) for a and —a. 

Figure 8: Overshoot for (J, h, a) in Regime (la) or for (J, /i, —a) in Regime (Id). Parameters: 
(J,/i) = (0.95,0.01), a = 0.46. 
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Figure 9: Absence of overshoot for (J, /i, a) in Regime (lb) or for (J, /i, —a) in Regime (Ic). 
Parameters: {J,h) = (0.3,0.04), a = 0.28. 
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A(m) 




Figure 10: Plot of A and intersection with the constant tanh(2/i). 

exists a unique = z~^{J,h) > such that k'^'^[z^) = and, in addition, if tanh(2^) < 
max A{m), then tliere exist —1 < < z ^ < such that k'^'^{z~) = k'^'^{z '") = (see 

me[-l,0] 

Fig. 0). 

Step B: Existence of m^, and relation between moo CLnd a. 

(Ba) Existence of m'^ : k'^'^{a) > together with a > and Step A imply a < z^. Since 
<I>(a) > and <I>(z+) < 0, it follows that there exists a such that < a < < z^ and 



= 0. 



(3.8) 



The latter in turn implies that k'-''^ (m'j^)'^ = — a^. This, together with k'^'^{m'^) > 0, 
implies that k'^'^{m~^) < m^, which leads to m°° < m^. 

(Bd) Existence of m^: As in (Ba), ^{z~) < and <l>(a) > imply that there exists a 
such that z" < < q < and 

^.(m^) = 0. (3.9) 



(Be): k'^'^{a) > and a < imply a < moo- This follows from the fact that k'^'^{a) > a 
implies a < rrioo by (|l.33| ). 

(Bh): k-''^{a) < and a > imply a > nioo- Again, this is a consequence of ( |1.33| ). 



Step C: Consequence of the positivity of times. Only positive solutions of Equation ( |3.4| ) are 
of interest. This implies the constraints 



r]F{m) := mk-^'^{m) + {a'^ -m'^) + \a\y^ 



tpim) > <?= 
and 

tUm) > ^ 
The functions rjp and r]L satisfy 

r]F{a) = ak'^'^{a) + \a\\k'^'^{a 



— j>0 if m2>a2, 
< if m"^ < a^. 



/ N , jh . s , , , / — 7 — r I > if > a^, , 

r]L(m) := m k'''^im) + (a^ - rr?) - a J^im) { ^ ' (3.11 

< if m"' < a^. 



'>0 if a F'''(a) > 0, 
= if a A:^''*(a) < 0, 



(3.12) 
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r^L{a) = ak-''^{a)-\a\\k''\a)\ = { ' (3.13) 

I < if a ^'^{a) < 0. 



Also, from ( pTssl) , 



0, 



0. 



(3.14) 



Last line implies that is a root of t^l but there is no change of sign aromid it. Finally, 
from expressions ( |3.1C1| )-( 3.14 ) we conclude that: 



The zeros of rjp, rji are a subset of {m°°, ita}. 

The intervals in which r]p and ry^ satisfies the constrains ( |3.10| )-( 3.11 ) are: 



Regime 


(la) 


(16) 


(Ic) 


(Id) 


Condition 


a > 0, y^^{a) > 


a > 0, A;'^'''(a) < 


a < 0, A;-^'''(q) > 


a < 0, k-^^^{a) < 


> 


[a, mj] 










> 




, a] 


[a, m,°°] 





We observe that in regime (la) each value of m E [a Amao,Tn'^j is attained at two different 
times tp ("First") and tp ("Last"). The same happens in regime (Id) for m G [m/?— , aVrrioo]- 
These phenomena correspond respectively, to an over and an under shoot. The proof is 
completed by showing that the trajectories have the right monotonicity properties. 
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Figure 11: r]F (dotted) and rji (dashed) for (J, /i, a) in regime 1(a) and for (J, h, —a) in regime 
1(d). Parameters: ( J, /i) = (0.95,0.01), a = 0.46 as in Fig. |. 



Step D: Monotonicity. 

By using implicit derivation we get 



d 



m 



dt 



dh,a 



m 



2 csch(2t)|mcsch(2t) - acoth(2t)| 
[kJ'^']'{m) - coth(2t) 



(3.15) 



On the other hand, if m = m{t) is a critical point (fc'^''^(rfi) = /(^Q,(m)), then <I>(m) 
(rncsch(2t) — acoth(2t))^. Hence, 
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dm , . 



= 0. 



Splitting in cases according to different values of m, it is not hard to conclude from ( |3.8| ) and 
(|3.15| ) the following monotonicity properties of the trajectories. 



Regime 


(la) 


(lb) 


(Ic) 


(Id) 


rh incr. 


< t < t(m+) 





< t < oo 


t{m^) < t < oo 


rh deer. 


t{ni^ < t < oo 


< t < oo 





<t < t{m~^) 



This concludes the proof of part 1 of Theorem 1.9 



3.2 Bifurcation 

Bifurcation proofs rely on the following facts. 

(Bl) For short times there is a unique critical point, close to a. 

(B2) Therefore, in order for bifurcation to occur a local maximum and a local minimum must 



appear in the course of time. Given condition ( 1.22 ), usual arguments imply that two 
(or more) stationary points appear at times larger than t if the curves and k'^'^ 
become tangent at a certain magnetization rh. The pairs (m, t) are determined by the 
following two equations (a similar argument was used in the proof of Theorem |1.7K iii)): 



[A;-^''^]'(m) =coth(2t), 

k'''^{7h) = mcoth(2f) - acsch(2t). 

Inserting the first equation into the second, we get 

m\k'''^^^\fh) — k''''^{rh) 



(3.16) 



F{m) :-- 



csch [arccoth( [fc"''''] '(m))] 



a. 



(3.17) 



We are left with the task of determining whether or not this equation has solutions. 
Note that 



F'{m) = 



[k-^'^] "(m) =0 or m = ^'^{m) [y^^] '(m) 



(3.18) 



In what follows all the assertions about F can be checked by using the equivalence in 
( 3.1^ ) and doing a straightforward analysis of k^'^. 



(B3) t I— 7- Ct,a is continuous with respect to || • ||oo- Hence, when a new minimum appears it 
cannot be a global one. Both a i— )• Ct^a and h i— )• Ct^a are also continuous. 

(B4) When t — )• oo we have two global minima ibm°° if /i = 0. If /i > the symmetry is 
broken and there is only one global minimum m°° > 0. 



25 



Whenever a local maximum/local minimum appears (disappears), we will refer to this 
behaviour as LMLMA (LMLMD). We proceed by looking at /i = and h ^ separately. 

Part (2) (h = 0, see Fig. O). The scenario for a = has already been proven in Theorem 



1.7. We concentrate on a > 0; this is no loss of generality due to the antisymmetry of F. 



Claim: Whenever a > 0, negative solutions of ( 3.17 ) can not cause bifurcations. In fact 



let tnc be the time at which the critical point in the negative side emerge. Let also, 

d{t) := Ct,a(m_(t, a)) - Ct,a{rh{t, a)), t> U, 



where m_(t, a) is the negative local [because of (B3)] minimum and m{t,a) is the global 
minimum of Ct^a- The last one is positive due to (Bl). and the supposition of a > 0. By 
definition, d{tnc) > and by (B4) liuit^oo d{t) = 0. Doing calculations similar to ( |2.15D and 
using that m_(t,a), m{t,a) are both critical points, we get that d'{t) < V t > t„c. This 
proves the claim. 



In what follows we focus in equation (p. 17 ). Owing to the previous claim, in order to 
bifurcation to occur a positive solution of (|3.17) is needed. 



(2a-b) If < J < |, then F'{m) < for all m G (-1,+1). Hence, for all a > there is 



only one solution of ( 3.171) . This solution turns out to be negative and hence it can not 



correspond to a bifurcation. 

(2c) If J > |, F has only one global maximum on the positive side, with value Ub = 
Ub{J) > 0. Combining B1.-B4., we get that there is bifurcation if and only if a S 
[0,t/B]=Im(F|[o,i]). 



Part (3) {h > 0, see Fig. 
Remark: If /i > 0, then B4. (Symmetry breaking) allows the appearance of solutions of (|3.17] ) 
leading to bifurcations. 

Once more, let study the different scenarios for F when h > 0. 

(3a) If < J < 1, then Im(F) C [-1, +1]^ Therefore (|3l7| ) has no solution for any \a\ < 1 
and there is no bifurcation. 

(3b) If 1 < J < |, then F has a unique maximum for m G [0, 1] with value Ub = Ub{J, h) < 0, 
and [— 1,C/b] = Im(-^|[o,i])- Arguing as in the claim of Part 2 for a > Ub, we conclude 
that there is bifurcation if and only if a G [—1, Ub]- 



(3c) Assume J > 



2- 



For h > small enough the behaviour is "close" to the h = case due to the 
continuity of /i i— t- Ct^a with respect to the infinite norm. 

Indeed, there exists Lb := min[_i g] ~ —Ub{J,0) and Ub ■= max[o^i] F ^ 
Ub{J,0), with {—1,Ub] = Im(F|[o,i]). There are different regimes for a: 

(a) For a < Lb, there is a unique solution of ( |3.17 ), which is in the positive side, 
leading to a bifurcation [because of (B4) ] . 

(b) For > a ^ Lb there are both a negative and a positive solution to ( |3.17] ) . 
Both lead to bifurcations, the negative one by continuity B3. and the positive 
due to B3.-B4.. The negative solution appears earlier in time. We write 
SB = SB{a) for the time of the first (negative side) bifurcation and tB = is (a) 
for the second bifurcation time. 
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Figure 12: First row: Plot of m i— F{m). % = LMLMA not leading to bifurcation; ▲ = 
LMLMD; ■ = LMLMA leading to bifurcation; = bifurcation . Second row: Plot of 
m Ct^aim) for different times. Short time (dotted) bifurcation (solid), long time (dashed). 



(c) For < a ^ Ub there is LMLMA on the positive and negative sides. As in 
the /i = case, the negative one does not lead to a bifurcation, and thus only 
one bifurcation occurs, which happens to be in the positive side. 

(d) By the continuity property (B3) and the monotonicity of a i— t- sb{oi) (proved 
below) , the two previous regimes coalesce, leading to an intermediate value 
Mt £ {Lb, Ub) such that trifurcation occurs at a = My. 

(e) For a > Ub there is no positive solution to ( 3.1'<j ) with a G [0,1]. Hence, no 
bifurcation occurs. 



2. The limit /i oo in ( |3.17] ) yields 



^a'^'(m) + b-^'{m) 



a. 



Hence, for h > large enough we get a behaviour similar to (3b), but with 
UB{J,h) > 0. 

3. The existence of h* follows from the continuity of the function h i— )• Ct^a commented 
in (B3) with respect to h. 
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(3b), /i small, J = 1.42, /i = 0.15 (3b), /i large, J = 1.42, /i = 1.6 




(3c), h small, J = 2.9, h = 0.15 (3c), /i large, J = 2.9, /i = 1.6 



Figure 13: Plot of m i— F{m) for different regimes of J when h > 0. # = LMLMA not 
causing a bifurcation; A = LMLMD; ■ = LMLMA causing a bifurcation. 



3.2.1 Monotonicity of the functions is (a) and SB{a) 

The bifurcation times are characterized by the following equations: 

y'^{m2) = ltsA^2), (3.19) 
CtB,a{rhi) = Ctg,ai'm2)- 

The first two equations say that rhi and m2 are stationary points at the same time tB, while 
the third one establishes the equality of costs at this time ts- Taking the derivative with 
respect to a of the third equation we get 

_ da da ^3 20) 

A straightforward computation using the first two equations shows that < 0, which implies 
that a I—)- tsia) is continuous and decreasing. A similar argument shows that a 1— )• SB{a) is 
continuous and increasing. 
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